3D XY and lowest Landau level fluctuations in deoxygenated YBa 2 Cu 3 07_j thin films 
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Conductivity measurements reflect vortex solid melting in YBa2Cu307-<s films. Field- 
independent glass exponents v g ~ f .9 and z g ~ 4.0 describe the transition T g (H) for < H < 26 T. 
At low fields, 3D XY exponents vxy — 0.63 and zxy — 1-25 are also observed, with zxy smaller 
than expected. These compete with glass scaling according to multicritical theory. A predicted 
power-law form of T g {H) is observed for 0.5T C < T g < T c . For T g < 0.5T C , 3D XY scaling fails, but 
a lowest Landau level analysis becomes possible, obtaining T C 2(H) with positive curvature. 
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The nature of fluctuations near the superconducting 
to normal state transition in high-temperature supercon- 
ductors (HTSCs) is still a matter of controversy. Sev- 
eral distinct fluctuation types and regions have been pro- 
posed, e.g., 3D XY fluctuations at low fields, lowest Lan- 
dau level (LLL) fluctuations at high fields, and glass- 
like fluctuations (for disordered HTSCs) near the finite- 
field transition T g (H). However, experimental analyses 
based upon the different scaling theories lead to conflict- 
ing results. This situation is most evident for competing 
3D XY and LLL fluctuations, both of which are sup- 
ported experimentally, in the same region of the phase 
diagram, in spite of being incompatible |l[] . 

The 3D XY transition is driven by phase fluctuations 
of a complex order parameter (OP), which fall into the 
universality class of the A transition in 4 He. The zero- 
field, "intermediate" (nonelectrodynamic) phase fluctu- 
ations of the HTSCs are thought to be of this type 
[§. At T = T c (and H = 0), these fluctuations di- 
verge in size, driving the resistive phase transition. Re- 
cent experimental evidence supporting this picture is 
found in specific heat magnetization Q^], penetra- 
tion depth ||, and current-voltage (J-E) measurements 
The finite-field transition T g (H), which is similarly 
driven by phase fluctuations of the OP, joins smoothly 
to T c = T g (H = 0). However, the glass and 3D XY fluc- 
tuations exhibit distinct scaling functions and exponents 

Fluctuations of the OP amplitude occur near the upper 
critical (mean-field) temperature T C 2(H). These fluctua- 
tions drive the Cooper pair density to zero, but do not 
correspond to a true transition; superconducting order 
vanishes at the slightly lower temperature T g (H). In 
the low-field region, the distinction between OP ampli- 
tude and phase fluctuations results in the dominance of 
3D XY critical behavior. At high fields, this distinction 
is not present, yielding a different type of behavior, most 
conveniently described in terms of the Ginzburg-Landau 
LLL approximation, with its corresponding scaling the- 
ory ^7^| ■ Experimental evidence in support of this be- 
havior is found in specific heat p]Tt2j, magnetization 



|10| , |12| , and J-E characteristics | [10jJl3j . A crossover is 
expected between the low-field (3D XY) and high-field 
(LLL) behaviors, and its clarification is fundamental in 
the investigation of HTSC fluctuations [Q . 

It may appear that glass fluctuations near T g (H) only 
complicate the story, since they compete with both the 
low- and high-field fluctuations. However, in this letter 
we suggest, to the contrary, that glass fluctuations help 
to identify low- and high- field behaviors, through the use 
of multicritical scaling theory. Such theories are appli- 
cable when fluctuations of different types compete for 
dominance. For example, in the low-field limit, T g (H) 
joins T c = T g (H = 0), forcing the distinct fluctuation 
types to coexist near the multicritical point T c Q. In 
this paper, multicritical predictions augment the 3DXY 
and LLL theories, thus clarifying their applicability to 
HTSCs. 

To address these issues, it is desirable to work in both 
the low- and high- field regions of the phase diagram. In 
YBa 2 Cu 3 7 -5 samples, H c2 (0) is very large (> 100 T) 
for "optimal" (5 ~ 0.05) stoichiometry. We therefore 
focus on deoxygenated films, for which magnetic field 
scales are relatively small, allowing H C 2(T) to be accessed 
over a wide temperature range. Three "optimal" c-axis 
YBa2Cu307_,5 films, approximately 4000 A thick, were 
prepared, and then deoxygenated, as described elsewhere 
p"5| . Films with resulting stoichiometries of 5 ~ 0.24, 
0.57, and 0.59 (±0.05) were produced, corresponding to 
T c s of 77, 61, and 56 K, respectively. The films were 
patterned into 100 /.im x 2000 /.im bridges using laser 
ablation. Isothermal J-E and conductivity {a) curves 
were obtained using a conventional four-point geometry. 
As typical for underdoped samples, the normal state con- 
tribution could not be eliminated from a, as it is not yet 
well characterized. However, this background contribu- 
tion should not affect the present results greatly, in the 
temperature range of interest, due to the insulating na- 
ture of the normal state. Magnetic fields were applied 
perpendicular to the film surface (parallel to the c-axis) 
in the range < H < 26 T. 

Transition temperatures T g , corresponding to the con- 
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tinuous vortex solid melting transition, were deduced 
at each applied field using the scaling ansatz of Fisher, 
Fisher, and Huse (FFH) |2j. (Note that this ansatz is 
isomorphic with any two-exponent scaling theory of the 
melting transition.) The appropriately scaled conductiv- 
ity is given by (J/E)\T - T^^-i), while the scaled 
current is (J/T)\T — T g \~ 2lJg . Here, v g and z g are the 
static and dynamic glass scaling exponents, respectively. 
The H g {T) phase boundaries determined in this way are 
shown in Fig. 1. The obtained scaling exponents are in- 
dependent of 8 pl| and H, with values of v g = 1.8-1.95 
and z g = 4.0-4.1. This is consistent with the notion of a 
single, 3D, glass universality class. 

As a first application of multicritical scaling, we note 
that the form of H g (T) [the inverse of TJH)] is speci- 
fied at low fields by the 3D XY theory ||@: H g (T) = 
H*(l - T/T c ) 2vxY , where v XY v g ) is the 3D XY 
static exponent. The phase boundaries satisfy this rela- 
tion over a wide temperature range, 0.52^ < T < T Cl as 
shown in Fig. 1 (inset), identifying the crossover temper- 
ature T b ~ 0.5T C as the limit of 3D XY scaling. The 
slope of these curves gives a (sample-averaged) exponent 
vxy = 0.63 ± 0.04, which compares favorably with the 
expected value of 0.669 |ll|. This is then used to de- 
termine the sample-dependent field scale H* . For the 
films used here, H* is in the range 7-19 T, as compared 
to "optimal" samples, for which H* is on the order of 
50 T. The crossover field H b = H g (T = T b ) was stud- 
ied previously, and was suggested to separate 3D from 
2D behavior along H g (T) 0Q. However, the field- 
independence of the glass transition does not corroborate 
this conclusion. Below, we demonstrate instead that for 
fields H < H b , the 3D XY description is in good agree- 
ment with the data, while for H > H b , a self-consistent 
LLL description becomes possible. 

It is possible to determine T c at zero field, using 
the FFH ansatz However, this procedure is known 
to become uncontrolled at very low fields, obtaining 
surprising results, such as nonuniversal exponents |jl9| . 
Therefore, we develop here a more reliable "crossing 
point" scaling technique, by extending the 3D XY anal- 
ysis to finite fields P0|| . The subsequent field scal- 
ing hypothesis involves the 3D XY scaling variable 
|,|: x = (H*/H)V 2vxY (T - T c )/T c . The scaling 
of the ohmic conductivity an can then be written as 
a n (H/H*) ( - ZXY - 1 ^ 2 = s(x), for which the asymptotic 
behavior is known |Jj2(|: s(x) ^ (x + l)^^ -1 ' as 
x — > —1, corresponding to T —> T g (H); and s(x) ~ 
j.-"xy(«xv-i) asu +oo, corresponding to H — > 0. 

The crossing point method proceeds from the defini- 
tion of x: if H > and T — T c , then x — 0, and 
must therefore be independent of H . It follows that in 
this limit, a n (H / H*)( ZXY ~ 1)/2 should also be indepen- 
dent of H. Data sets an(T), obtained at constant fields 
H, must then all cross at T — T c when plotted as in 
Fig. 2, provided that zxy has been chosen correctly. As 



observed in the insets, the crossing point method places 
strong constraints on the exponent and transition tem- 
perature, which we identify as zxy — 1-25 ± 0.05 and 
T c = 60.8 ± 0.4, for the film shown. These results are 
corroborated by low field data (ff < 0.1 T) @. Since 
3D XY fluctuations are most prevalent near T ~ T Cl it is 
helpful to think of this method as optimizing scaling near 
this temperature. We emphasize that zxy obtained in 
this way disagrees with the expected Q diffusive dynam- 
ics (zxy — 2), and also with other experimental obser- 
vations p9| , pT[ |. It is our opinion that the crossing point 
method can achieve better results than those analyses, 
due to the incorporation of finite field data in the scaling 
procedure. 

The estimates for zxy-, ^xy, and T c , determined 
above, can be checked through a full scaling analysis of 
an, as shown in Fig. 3. For that film, the initial estimates 
could not be improved upon. The fitting is excellent at all 
low fields, becoming optimal at T — T c . For comparison, 
the best fit, using the expected exponents zxy — 2 and 
vxy = 0.669, is shown in the inset. As found in previous 
analyses scaling using zxy — 2 suffers at the lowest 
fields, and perhaps more importantly, near T c . The scal- 
ing results found here are therefore an improvement over 
previous analyses. 

The crossover between low- and high-field behavior, 
observed above using multicritical scaling, can be for- 
mulated more conveniently as follows. Since all sample 
dependence of the 3D XY scaling variable x is absorbed 
into the characteristic field H*, the divergence of an in 
Fig. 3, must occur at a universal value of x = — 1 [ p0[ . 
Multicritical self-consistency therefore requires the scal- 
ing variable evaluated along the phase bound- 
ary T g (H)] to remain field-independent in the 3D XY 
scaling region, as shown in the lower half of Fig. 4. Er- 
rors in the determination of x g are magnified at lowest 
fields, where the difference T c — T g (H) is small. Devia- 
tion from 3D XY multicritical self-consistency becomes 
apparent for H > Hi,. 

Several fluctuation models are candidates for describ- 
ing the upturn of H g (T) when H g > H b in Fig. 1 
(inset). Here, we consider the 3D LLL model, using 
multicritical theory to place restrictions on the allow- 
able scaling. The scaling technique is constructed in 
analogy with the preceding 3D XY analysis. In the 
LLL theory, a natural scaling parameter emerges 
y = {H*T C /HT) 2 / 3 [T - T c2 (H)}/T c , where T c and H* 
have been used here to make y dimensionless. In this 
analysis it is T C 2(H) which must be determined by scal- 
ing. Although T C 2{H) is often assumed to be linear 
this restriction becomes too severe for the underdoped 
samples used here. Instead, a LLL crossing point method 
is now constructed, which allows the form of T C 2(H) to 
be ascertained. As described below, this analysis ob- 
tains T C 2(H) curves with positive curvature — an inter- 
esting feature which has previously been associated only 
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with magnetically doped HTSCs 

The LLL crossing point method is described as follows. 
We make use of the 3D LLL scaling ansatz M , which we 
write as a n (HT? / iJ*T 2 ) 1/3 = F 3B (y). In analogy with 
the 3D XY case, an{HT^ /£TT 2 ) 1/3 must be indepen- 
dent of H when T = T c2 {H). To begin, a value of T c2 
is first assumed for a particular reference field Hq. The 
temperatures T C 2{H) consistent with this choice are then 
obtained for other fields (H ^ Ho). This is accomplished 
by plotting a n (HT? / H*T 2 ) 1/3 vs. T — T c2 (H) for (fixed 
H) <Jn(T) data sets, then adjusting T c2 (H) for each field, 
until a crossing occurs at T — T c2 (H) = 0, similar to 
Fig. 2. For fields H ^ H , T c2 (H) is thus a function of 
the original choice of T c2 (Ho), reducing the following fit 
to a single parameter. All T c2 {H) curves found in this 
way exhibit positive curvature. It is once again helpful 
to view the crossing point as a method for optimizing 
scaling in the most essential temperature region; in the 
LLL case this is near T c2 (H). 

We are now left with one fitting parameter, T c2 (Hq), 
which cannot be estimated from the crossing point 
method alone. The full LLL scaling procedure is now 
used to determine the single fitting parameter, while si- 
multaneously requiring multicritical self-consistency. In 
analogy with the 3D XY case, this means that the scaling 
variable y g [i.e., y evaluated at T — T g {H)] must remain 
field-independent in the LLL scaling region. We find 
that by using T c2 (H) obtained from the crossing point 
method, multicritical self-consistency cannot be met at 
low fields. Since the LLL theory has greatest justifica- 
tion at high fields, we attempt, instead, to obtain self- 
consistency in the high-field region. The outcome of the 
final scaling procedure is shown in Fig. 5. The (approx- 
imate) field-independence of y g is evident for the entire 
high-field range H > Hb, as shown in the top half of 
Fig. 4. T c2 {H) is obtained with only small uncertainty, 
as shown in Fig. 5 (inset). 

We comment finally on the difference between the 
present results and those of Refs. (To| and |fl3| . In 
our work, LLL scaling is found to be multicritically 
self-consistent only at high fields (H > Hi,), while 
Refs. jl(],[n|, which do not check for self-consistency, find 
that LLL scaling is applicable at low fields (H < Hb). 
(Note that Hb is very large in the "optimal" samples used 
by those authors.) We speculate that scaling could be ac- 
complished in Refs. (h]]l| only by allowing diminished 
scaling quality in precisely the region where the quality 
should be highest [i.e., near T c2 (H)]. In the present 
work, this situation is avoided by optimizing scaling near 
T c2 (H) from the outset. 

After the completion of this work, we learned of simi- 
lar 3D XY results, obtaining zxy ^ 2 from conductivity 
scaling [ p3| . 

We thank S. Girvin, A. MacDonald, S. Pierson, 
Z. Tesanovic, and especially P. Muzikar for many helpful 
discussions. This work was supported through the Mid- 



west Superconductivity Consortium (MISCON) DOE 
grant # DE-FG02-90ER45427, the Materials Research 
Science and Engineering Center (MRSEC) Program of 
the NSF under Award # DMR-9400415, and the Na- 
tional High Magnetic Field Laboratory (Tallahassee). 



[1] I. D. Lawrie, Phys. Rev. B 50, 9456 (1994). 

[2] M. P. A. Fisher, Phys. Rev. Lett. 62, 1415 (1989); D. S. 

Fisher, et al, Phys. Rev. B 43, 130 (1991). 
[3] S. E. Inderhees, et al, Phys. Rev. Lett. 66, 232 (1994); 

G. Mozurkewich, et al, Phys. Rev. B 46, 11914 (1992); 

N. Overend, et al, Phys. Rev. Lett. 72, 3238 (1994); M. 

A. Howson, et al, Phys. Rev. Lett. 74, 1888 (1995); N. 
Overend, et al, Phys. Rev. Lett. 75, 1870 (1995). 

[4] M. B. Salamon, et al, Phys. Rev. B 47, 5520 (1993); M. 

B. Salamon, et al, Physica A 200, 365 (1993). 

[5] M. A. Hubbard, et al, Physica C 259, 309 (1996). 

[6] S. Kamal, et al, Phys. Rev. Lett. 73, 1845 (1994); S. M. 

Anlage, et al, Phys. Rev. B 53, 2792 (1996). Note also 

the counterexample of Z.-H. Lin, et al, Europhys. Lett. 

32, 573 (1995). 
[7] N.-C. Yeh, et al, Phys. Rev. B 47, 6146 (1993); M. A. 

Howson, et al, Phys. Rev. B 51, 11984 (1995). 
[8] A. J. Bray, Phys. Rev. B 9, 4752 (1974); D. J. Thouless, 

Phys. Rev. Lett. 34, 946 (1975); G. J. Ruggeri and D. J. 

Thouless, J. Phys. F 6, 2063 (1976); R. Ikeda, et al, J. 

Phys. Soc. Jpn 58, 1377 (1989); S. Hikami and A. Fujita, 

Phys. Rev. B 51, 6379 (1990). 
[9] S. Ullah and A. T. Dorsey, Phys. Rev. Lett. 65, 2066 

(1990); Phys. Rev. B 44, 262 (1991). 
[10] U. Welp, et al, Phys. Rev. Lett. 67, 3180 (1991). 
[11] B. Zhou, et al, Phys. Rev. B 47, 11631 (1993); A. Junod, 

et al, Physica C 211, 304 (1993); S. W. Pierson, et al, 

Phys. Rev. Lett. 74, 1887 (1995); M. Roulin, et al, Phys. 

Rev. Lett. 75, 1869 (1995); M. Roulin, et al, Physica C 

244, 225 (1995); S. W. Pierson, et al, Phys. Rev. B 53, 

8638 (1996). 

[12] O. Jeandupeux, et al, Phys. Rev. B 53, 12475 (1996). 

[13] D. H. Kim, et al, Phys. Rev. B 45, 10801 (1992). 

[14] Z. Tesanovic, Phys. Rev. B 51, 16204 (1995). 

[15] L. Hou, et al, Phys. Rev. B 50, 7226 (1994). 

[16] J. C. LeGuillou and J. Zinn- Justin, J. Phys. (Paris) 46, 

L137 (1985). 
[17] J. Deak, et al, Phys. Rev. B 51, 705 (1995). 
[18] A. Schilling, et al, Phys. Rev. Lett. 71, 1899 (1993); C. 

C. Almasan and B. Maple Phys. Rev. B 53, 2882 (1996); 
K. Kishio, in Physical Properties of High Temperature 
Superconductors, edited by D. M. Ginsberg (World Sci- 
entific Singapore, 1996). 

[19] J. M. Roberts, et al, Phys. Rev. B 49, 6890 (1994). 
[20] M. Friesen, et al, unpublished. 

[21] J. C. Booth, et al, Phys. Rev. Lett. 77, 4438 (1996). 
[22] Yu. N. Ovchinnikov and V. Z. Kresin, Phys. Rev. B 54, 
1251 (1996), and references therein. 



3 



[23] D. Ginsberg, J.-T. Kim, and N. Goldenfeld (unpub- 
lished) . 




T(K) 

FIG. 1. Superconducting phase boundaries for three de- 
oxygenated YBa2Cu307_,5 thin films. Lines are a guide 
to the eye. Inset shows the same phase boundaries plot- 
ted logarithmically. Power law behavior is evident for 
H <H b . 
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FIG. 2. 3D XY crossing point method, shown for the 
8 = 0.57 film. Symbols correspond to different fields: 1, 
2, 3, 4, 5 T. The appropriate choice of zxy causes lines 
to cross at a single point, identifying both zxy and T c . 
Poor crossing behavior is observed for (slightly) incorrect 
values of zxy (insets). 




-1 -0.5 0.5 1 1.5 
(H*/H) 1/2v « (T-TJ/T. 



FIG. 3. Full 3D XY scaling of the ohmic conductivity 
an for the 8 ~ 0.57 film. Scaling is successful only for low 
fields (H < Hb ^ 5.1 T). Expected asymptotic behaviors 
are shown as solid lines. Insets show the poor scaling 
obtained using the expected exponents vxy = 0.669 and 
z X Y = 2. 
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FIG. 4. Multicritical self-consistency: 3DXY and LLL 
scaling parameters x g and y g , respectively, evaluated at 
T = T g (H). Dashed line indicates the crossover Hi,. Scal- 
ing variables should remain field-independent in their re- 
spective scaling regions. For H < H b , 3DXY scaling is 
self-consistent, while for H > Hb, LLL scaling is self- 
consistent. 
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FIG. 5. Full LLL scaling of cr n for the 8 ~ 0.57 film. 
Multicritically self-consistent scaling is successful only 
for high fields (H > H b ~ 5.1 T). Inset shows T c2 {H) 
and Tg(H) for the same film. [T c2 (H) is speculative for 
H < H b .} 
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